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ABSTRACT 

We present a new, publicly available model of the extragalactic background light (EBL) and corresponding 
77 opacity for intergalactic gamma-ray absorption from z = up to z = 10, based on a semi-analytical model of 
hierarchical galaxy formation that reproduces key observed properties of galaxies at various redshifts. Including 
the potential contribution from Population III stars in a simplified way, the model is also broadly consistent with 
available data concerning cosmic reionization, particularly the Thomson scattering optical depth constraints 
from WMAP. In comparison with previous EBL studies up to z ^ 3-5, our predicted 77 opacity is in general 
agreement for observed gamma-ray energy below 400/(1 +z) GeV, whereas it is a factor of ^ 2 lower above 
this energy because of a correspondingly lower cosmic star formation rate, even though the observed UV 
luminosity is well reproduced by virtue of our improved treatment of dust obscuration and direct estimation 
of star formation rate. The horizon energy at which the gamma-ray opacity is unity does not evolve strongly 
beyond z ^ 4 and approaches ^ 20 GeV. The contribution of Population III stars is a minor fraction of the EBL 
at z = 0, and is also difficult to distinguish through gamma-ray absorption in high-z objects, even at the highest 
levels allowed by the WMAP constraints. Nevertheless, the attenuation due to Population II stars should be 
observable in high-z gamma-ray sources by telescopes such as Fermi or CTA and provide a valuable probe of 
the evolving EBL in the rest-frame UV. 

Subject headings: cosmology: diffuse radiation - gamma rays : theory 



1. INTRODUCTION 

The extragalactic background light (EBL), the diffuse, 
isotropic background radiation from far-infrared (FIR) to ul- 
traviolet (UV) wavelengths, is believed to be predominantly 
composed of the light from stars a nd dust integrated oyer the 
entire history of the Universe (see iDwek & KrennricEl2012L 
for reviews). The observed spectrum of the local EBL at z = 
has two peaks of comparable energy density. The first peak 
in the optical to the near-infrared (NIR) is attributed to di- 
rect starlight, while the second peak in the FIR is attributed to 
emission from dust that absorbs and reprocesses the starlight. 

The precise determination of the EBL has been a difficult 
task. Direct measurements of the EBL in the optical and NIR 
bands have been hampered by bright foreground emission 
caused by interplanetary dust, the so-called zodiacal light (see 
iHauser & D wek 200 1 , for reviews). Recently, Matsuoka et al. 
(1201 ll) reported measurements of the EBL at 0.44 fim and 
0.65 fj,m from outside the zodiacal region using observational 
data from Pioneer 10/11. On the other hand, integration over 
galaxy number counts provide a firm lower bound on the EBL, 
and the observed trend of the counts with magnitude indi- 
cates that the EBL at z = has been largel y resolved into dis- 
crete sources in th e optical/NIR bands (e.g.fMadau & Pozzettil 
HOOO; Tota nTetani200llKeenan et alJl2010^ ■ even when the 
effect of incomplete detection due to cosmological dimming 
of surface brightness is taken into account ( Totani et al. 2001). 

The EBL can also be probed indirectly through observa- 

|E-mail: yinoue@slac.stanford.edu] 



tion s of high-en ergy gamma rays from extragalactic objects 
(e.g. lMazin & Ra ue 2007). Gamma-rays propagating through 
intergalactic space can be attenuated by photon-photon pair 
production interactions (77 — > e^e~) with low-energy photons 
of the EBL. For gamma-rays of given energy E^, the pair pro- 
duction cross section peaks for low-energy photons with en- 
ergy 
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where is the electron mass and c is the speed of light. In 
terms of wavelength. Apeak = 1.24(£'^[TeV])/im. Measuring 
the resultant attenuation features in the spectra of extragalac- 
tic GeV- TeV sources would offer a valuable probe of the EBL 
that is indirect, yet unique in being redshift-dependent. Al- 
though this method can be limited by incomplete knowledge 
of the intrinsic spectra of the source before attenuation, by 
assuming a plausible range for such spectra, observations of 
blazars by current ground-based telescopes have been able to 
place relatively r obust upper limits to the EBL at z = an d 
up to z 0.5 (e.g. lAharonian et al.ll2006al: I Albert et"ani2008h . 
Combined with the lower limits from galaxy counts, the total 
EBL intensity at z = from 0.1 /im to 1000 fim is inferred 
to lie in the range 52-99 nW/m^/sr (Horiuchi et al. 2009). 
However, the EBL at higher redshifts is still highly uncertain. 

Currently available theoretical models for the EBL can be 
broadly categorized into three types. First, in backward evo- 
lution models, one starts from the observed properties of 
galaxies in the local Universe and describes their evolution 
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by extrapolating backwards in time in a par ameterized fash- 
ion (e.g. [M alkan & Stecker 1998; Totani & Takeuchil 120021: 
IStecker etalji2006i; f ranceschini et al.ii2008i) . This extrapo- 
lation entails uncertainties in the properties of the EBL that 
inevitably increase at high redshifts. Nevertheless, based on 
the observed, rest-frame K-band l uminosity functio r i (LF 
of galaxies from z = up to z = 4, iDommguez et al.l (l2011f 
we re able to model the EBL without any assumptions fo r the 
LF. lHelgason & Kashlinskvl(l2012l) : IStecker et al.l (120121) con- 
structed evolving EBL models in a relatively robust way by 
utilizing multiwavelength photometric survey data. 

Secondly, in forward evolution models, the basis is a de- 
scription for the cosmic star formation history (CSFH), over 
which models for the spectral energy distribution (SED) of 
the stellar population are convolved to obtain the evolving 
EBL (e.g. Kneiske et al. 2004; Finke et al. 2010). However, 
such models cannot follow the detailed evolution of key 
physical quantities such as the metallicity and dust content, 
which can significantly affect the spectral shape of the EBL. 
Furthermore, altho ugh most forward evolutio n models em- 
ploy the CSFH of iHopkins & BeacomI (120061) . it is known 
that this CSFH model overproduces the stellar mass den- 
sity ([Choi & Nagamine 2012), and is also in consistent with 
the ob served rate of core-c ollapse supernovae (jHoriuchi et alJ 
20111). Recent s tudies by 'Kobayashi et al.' (120 12|) show that 



Hopkins! ( 120041) ; IHopkins & Beacom (2006) may have over- 



estimated the CSFH at z > 1, arising from overcorrection for 
dust obscuration effects and in conversion from luminosity to 
star formation rate. 

Finally, rooted in the modern cosmological framework of 
large-scale structure formation driven by cold dark matter, 
semi-analytical models of hierarchical galaxy formation 
account for the merging history of dark matter halos as 
well as the physical evolution of the baryonic compo- 
nent, including the effects of gas cooling, star formation, 
metal enrichment, fe e dback h eating, etc. (Primack et al. 
2005t iGilmore et all l2009t lYounger & Hopkins .2011: 



Gilmore et al.l 12012k Such models successfully reproduce 



various observed properties of galaxies from the local Uni- 
verse up to z 6 (see e.g., Kauffmann et al. 1993; Cole et al. 
19941: jNagashima et al. 1999; Somerville & Primack 
iNagashi ma & Yoshii .2004: .Baugh et aL |2005| 
Nagashima et al 
Somerville et al. 
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120121) . At present, semi-analytical models 



can be considered the most detailed and well-developed 
models for the EBL over a wide range of redshifts. 

A subject that has yet to be fully explored in the context of 
the EBL and gamma-ray absorption is the epoch of cosmic 
reionization above z ^ 6. Measurements of NIR absorption 
troughs in the spectra of high-z quasars, together with those 
of anisotropics in the polarization of the cosmic microwave 
background (CMB), prove that the majority of intergalactic 
hydrogen in the Universe has been r eionized somewhere be - 
tween z ^ 1000 and z ^ 6 (see e.g. iRobertson et alJl20T^ . 
Although the most widely suspected source of reionization 
is UV photons emitted by early generations of massive stars, 
the observational constraints are still very limited, so that the 
actual sources, history and topology of cosmic reionization 
remain largely unknown. A closely related topic is the possi- 
ble existence and formation history of Population III (Pop- 
III) stars, very massive stars that are expected to originate 
in nearly metal-free conditions, particularly for the very first 
generation of stars appearing in the Universe, and th eir po- 
tential role in cosmic reionization (see e.g. .Fernandez et alj 



I20TI) . 

Since the current observational constraints on reioniza- 
tion all concern the neutral or ionized intergalactic gas, it 
would be very valuable and complementary to obtain inde- 
pendent information on the evolving, UV intergalactic ra- 
diation field itself. A unique and promising approach may 
be offered by the effects of gamma-ray absorption in-situ in 
high-energy sources at z > 6. UV radiation fields with suffi- 
cient intensities to cause cosmic reionization may induce sig- 
nificant gamma-ray absorption at observer energies above a 
few tens of GeV (Oh 2001; Inoue et al. 2010a). Based on 
a semi-analytical model of galaxy formation that includes 
Pop-Ill stars and reproduces a variety of r eionization-related 
ob servations ([Choudhurv & Ferraral 120061) . the recent study 
by llnoue et al.l ( 120 lOa ^) suggested that appreciable attenua- 
tion may be expected above ^12 GeV at z ^ 5 and down 
to ^6-8 GeV at z > 8-10, mainly caused by Pop-II stars 
at these epochs. However, the relative contribution of Pop-Ill 
stars was found to be difficult to discern observationally. On 
the other hand, without addressing the implications for reion- 
ization, some studies have concentrated on the prospects for 
constraining Pop-Ill star formation through gamma-ray ab- 
sorption in objects at lower redshifts (e.g. Gilmore 2012^. 

The Fermi gamma-ray space telescope (lAtwood e t al. 20(M 
Ferm i) may eventually detect blazars at z > 6 ( Inoue et aU 
I2OI Ih . and the Cherenkov Telescope Array dActis e t al. 201 J 
CTA) m ay possibly do the same for gamma-ray bursts (GRBs) 
(llnoue et al. in preparation). Therefore a deeper investiga- 
tion into the above issues is worthwhile and timely. The 
above studies (Oh 2001; Inoue et al. 20 10a; Gilmore |20 121) 
have not accounted consistently for the EBL resulting from 
ga laxy formation at lo wer redshifts. For example, the model 
of llnoue et all (12010 a) was optimized to describe the reion- 
ization epoch and did not include the contribution from Pop- 
I stars or dust, and thus could only evaluate the gamma-ray 
opacity above z = 4. 

In this paper, we present a new study of the EBL and con- 
sequent gamma-ray opacity, covering the entire redshift range 
from z = 0uptoz=10 within a consistent framework, ac- 
counting for the process of cosmic reionization, and includ- 
ing Pop-Ill stars in a simplified way. As a baseline model, 
we ado pt the Mitaka mod eQ of semi-analytical galaxy for- 
mation (INagashima & Yos hii 2004). The model can repro- 
duce various observed properties of galaxies such as the lu- 
minosity function (LF), lumino sity dens ity (LD), and stellar 
mass density (Nagashima & Yoshii 2004), as well as the LFs 
of high-redshift Lyman-break galaxi es (LBG s) and Lyman-a 
emitters (LAEs) up to z 6 (Kob ayashi et a l. 2007, 2010). 
As regards Pop-Ill stars, in view of the presently large theo- 
retical uncertainties on their formation efficiency, metal yield, 
conditions for transition to Pop-II star formation, etc., we 
choose not to fully incorporate them into our semi-analytical 
scheme. Instead, their formation history is characterized in 
a simple, parameterized way, which we constrain by model- 
ing the cosmic reionization process and comparing with avail- 
able observations, particularly the Thomson scattering opti- 
cal depth measured by the Wilkinson Microwave Anisotropy 
Probe (WMAP). 

We introduce our semi-analytical model in ^ Cosmic 
reionization is modeled and compared with observations in 
SI SI presents the results of our EBL models. The con- 

' Named after the city of Mitaka where the model was mainly developed 
at the National Astronomical Observatory of Japan. 
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Fig. 1. — Cosmic star formation liistory. Solid curve shows the to- 
tal in the baseline Mitaka model, while the dashed, dotted, dot-dashed 
and double-dot-dashed curves show the fractional contributions from stars 
with metallicity Z < 10^ = 5 X lO^'Z©, 5 X IO^'Zq < Z < 2 x IQ-^Zq, 
2 x 10"^ Zq < Z < Zq, and Zg < Z, respectively. We also plot the ob- 
servational dtoac(OTi^^ (2004), that deduced from LBGs 
lOuchi et al."2004l:IBou wens et al.,,20Q7,; Mannucci et al. 2007; Verma et aT] 
pSbV; Bouwens et al. 2008), and that inferred from GRBs I Yuksel et al. 200l; 
iKistler et al. 2009). For the LBG sample, limiting luminosities adopted by 
each author are indicated in the corresponding legend in term of the charac- 
teristic luminosity L, of the luminosity function. 



sequent gamma-ray opacity and comparison with current 
gamma-ray observations are described in ^ We conclude 
in ^ Throughout this paper, we adopt the standard cos- 
mological parameters of (h^riM.^A ) = (0.7,0.3,0. 7), and a 
Salpeter initial mass function (IMF: lSalpeteii[T955h within a 
mass range of 0.1 - 60 Mq. 

2. SEMI- ANALYTICAL GALAXY FORMATION MODEL 

In the framework of the Mitaka semi-analytical model of 
galaxy formation, we compute the evolution of baryon com- 
ponents within dark matter halos following the merger his- 
tory of dark matter halos. The time evolution of baryons 
in halos is computed with physically motivated, phenomeno- 
logical prescriptions for radiative cooling, star formation, su- 
pernova feedback, chemical enrichment, and galaxy mergers. 
Then, we can derive various physical and observational quan- 
tities for galaxies as well as the global average over the Uni- 
verse at any redshift, such as the CSFH, LFs, and dust con- 
tent of galaxies. A mock catalog of galaxies can be gener- 
ated, to be compared with vario us observations. Detai l s of th e 
Mitaka model are described in iNa^ashima & Yoshiil JIOOl : 
Eobayashi et al. ( 20071 IMol) . 

Since there are several free parameters in the phenomeno- 
logical prescriptions for baryons, we fix them to fit various 
observational data of the local galaxies such as B-band and 
K-band LFs, neutral gas fraction, and ga s mass-to-luminosity 
ratio as a function of B-band luminosity (INagashima & Yoshiil 
|2004). For consistency, we keep these parameters unchanged 
in this paper 

2.1. Cosmic Star Formation History 

Fig. [T] shows the CSFH expected in the Mitaka model 
at z = - 18. The CSFHs at various metallicity bins are 
also shown. Pop-Ill stars are expected to form from gas 
with metallicity below a critical value, such that the gas can 
only cool rather inefficiently through rotational transitions of 
molecular hydrogen, which leads to fragmentation into rela- 
tively massive protostellar clouds, and ultimately the forma- 
tion of very massive stars. Once the metalHcity exceeds this 



value, the gas can cool more efficiently via metal transitions, 
and a transition to the formation of less massive. Population 
II (Pop-II) s tars is thought to take place jMackev et al.ll2003l: 
iBromm & L oeb 2003; Yoshida et al. 20(^. However, the ex- 
act value of this critical metallicity has been hotly debated, 
ranging from Z = lO '^Zp = 10"^^ (ISchneider et aL.2006.) to 
Z = IO-^Zq = 10-3^ dAvkutalp & SpaansllMTH In this pa- 
per, we consider stars with metallicity Z < 10~* to corre- 
spond to Pop-Ill stars. We adopt a Salpeter IMF in the mass 
range of 0.1 -60Mq for all types of stars. Recent radiation- 
hydrodynamics simulations of Pop-Ill star formation suggest 
that their masses may be limited to < 4QMq due to radiative 
feedback effects (iHosokawa et al.ll201 li) . which would be in 
accord with our choice of the the maximum mass for Pop-Ill 
stars. 

We also plot the data compiled b y iHopkinsI (12004^ 
Pascale et al. (2009) : iRodighiero et all (2010); Karim et d] 
(i2011i): iCucciati et a L ('2012') , thaF deduced from LBGs 



Ouchi et alj l2004t Bouwens et all 120 07; Mannucc i et al.1 



120071; I Verma et al.H 2007; Bouwens etal. 2008, 2011b ), and 
that i nferred from GRBs (Yuksel etal. 2008; Kistler et al.l 
l2009h . Above z > 4, each LBG data points are obtained by 
integrating its LF down to a certain limiting luminosity which 
is parameterized by the characteristic luminosity L, of the LF. 
We note that different authors choose different limiting lumi- 
nosities. Our model shows the star formation in all galaxies, 
down to the faintest luminosities. 

There is an apparent discrepancy between our model and 
the observed CSFH data at 1 < z < 5. In our semi-analytical 
model, we directly estimate the star formation rate in a galaxy 
and evaluate the CSFH by integrating the star formation rates 
of all galaxies. Since the CSFH data are converted from the 
observed galaxy LFs, there are following observational uncer- 
tainties: faint-end slope of the LF, dust obscuration correction 
from UV data, contamination from old stellar populations to 
the IR luminosity, total IR luminosity modeling, and a conver- 
sion factor from luminosity to star formation rate. In the case 
of Hopkins (2004); Hopkins & Beacom (2006) CSFH model, 
there are inconsistencies with other observational informa- 
tion such as stellar mass density (e.g. Choi & Nagamine 2012h 
and co re collapse supernovae r ate iHoriuchi et al. 201 IJ). Re- 
cently, iKobayashi et al.l (120121) have shown that this discrep- 
ancy is due to the overcorrection of dust obscuration and star 
formation rate conversion. These uncertainties will result 
a fact or of ^2-3 overestimation of CSFH ( Kob avashi et al.l 
l2012h . Therefore, the discrepancy does not adversely affects 
our results, since our model can reproduce various observa- 
tions for galaxies. 

2.2. Stellar and Dust Emission 

We calculate the emissivity of the Universe at a given red- 
shift from the CSFH that is given as a function of metal- 
licity and the dust attenuation strength using stellar popu- 
lation synthesi s models . We u se stellar population synthe- 
sis models o f|Schaerer| (l200l for Z < lO""* (Pop-Ill) and 
iBruzual & C hariot ( 20031) forZ > 10"^ to calcu late the stellar 
emissivity. Although lBruzual & Charlo? (I2003h is commonly 
used and cover a wide range of metallic ity, the Z < lO"'* pop- 
ulation is not included. ISchaererl (l2003h can evaluate the stel- 
lar population evolution at Z < 10"^. For dust attenuation 



^ W e adopt Zp, ~ 0.02 ^ nders & Gre vessd [T989I : IGrevesse & Sauvall 
IT99I) . However. I^plund et alj f2009) updated the value of Z© ~ 0.0134. 
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Fig. 2. — Luminosity density at 1500 A. Solid line is the Mitaka model. 
The obser ved data at various redsh ifts are also shown as indicated in 
the figure <Wvder et alj|20 05: Schimi n]ovich et alj|200l IDahlen et aUi2007c 
IBouwens et aU2007l ; IRedd v et al. 200,^^ 



in galaxies and the inter galactic medium (I GM), ■we adopt 
iCalzetti et all (120001) and lYoshii & PetersonI (,19941 respec- 
tively. 

Various instruments have constrained the UV LD u p to z ^ 
6 ( Wvder et al. 2005; Schiminovich et al. 2005; Dahle n et all 
l2007t IBouwens et al.ll2007l: iReddv et al. 2008). Fig. |2]com- 
pares the UV LD from our model with observed LD data 
at a rest-frame of 1500 A. All the data agree well with our 
semi-analytical model. Comparison with other kinds of data 
are sh own in Nagashima & Yoshii (2004); Nagashim a et al.l 
(12005'); Kobayashi et al. (2007, 2010). 

Dust emits mid-IR (MIR) and FIR photons by reemitting 
the absorbed starh ght. For dust emission, we al so utilize the 
Mitaka model (see iMakiva et aLllin preparation! for details). 
We set the total IR emissivity is equivalent to the total ab- 
sorbed starlight energy by dust. Although our IR emissiv- 
ity model can rep roduce the local Herschel galaxy luminosity 
function (IVaccari et a l. 2010), the redshift evo lution of IR lu- 
minosity functi on (e.g. .Rodigh iero et al.''2010) has not been 
reproduced yet (iMakiva et aLl Tin preparation). Thus, we pre- 
dict an IR EBL at z = that is than limit of the EEL by a 
factor of two. To avoid this, we set our IR emissivity model 
three times more luminous than the usual version of the Mi- 
taka model. Our IR model prediction at MIR-FIR is uncer- 
tain, while our model is accurate at the UV-NIR band. We 
also note that gamma-ray opacity generated by dust emission 
is effective only above several tens TeV (see Eq. 

For the IR spectrum template, we use the IDale & Heloul 
(|2002) dust emission model. In their model, the IR spectral 
energy distribution (SE D) shape i s defined by the exponent 
a of Eq. 1 in .Dale & Heloul dlOOl . We adopt a = 1 .2 for all 
galaxies in our model to reproduce the peak wavelength of the 
FIR EBL. When we use a different a parameter, the position 
of the FIR peak wavelength will change. 

3. COSMIC REIONIZATION HISTORY 



Cosmic reionization is a major change in the ionization 
state of the Universe. After the recombination epoch, the Uni- 
verse was in the dark age where no radiation sources existed. 
As dark matter halos evolved from the first small inhomo- 
geneities in matter density field, some of baryons started to 
fall in the dark matter halo potential well. When the mass of 
the halo is high enough, the fallen gas formed stars and galax- 
ies. The first population of stars and galaxies can generate the 




Fig. 3. — Ionizing photon emissivity per comoving Mpc', dn\a^/dt, as 
a function of redshift. We set /esc = 0.2 for all models and observed 
data. Solid, dashed, and dotted curve shows the baseline, the lower-Pop-III, 
the upper-Pop-III model, respectively. Dot-dashed curve shows the ioniz- 
ing photon emissivity from Pop-Ill population in the baseline model. The 
fille(l symbol shows t he dat a derived from galaxy UV LFs I Yoshida et aQ 
12006; Bouwens et al? 12003; lOuchi et al. 2009). The data derived from 
the combination of hydrodynamical simulations and Lyof forest opacity 
( Bolton & Haehnelt 2007) are shown by open symbols. Thin dashed lines 
plot the model predictions of the ionizing photon emissivity that is required 
for maintaining hydrogen ionization in IGM ( Mada u et al- 1999.) with clump- 
ing factors of C = 1,3, and 10, from bottom to top. 



UV radiation. This rad iation would ionize the Universe (see 
iBarkana & Loebir200lL for reviews). In this section, we fol- 
low the reionization history of the Universe with our Mitaka 
model. 

3.1. Ionizing Photon Emission Rate 

The emission rate of hydrogen ionizing photon is required 
to follow the cosmic reionization history. Hereafter, we take 
the escape fraction of ionizing photons from a galaxy /esc at 
al l redshift to be 0.2 motivated by the numerical simulations 
of lYajima etaH (|2009l .2011). They also showed that /esc de- 
pends on halo mass. Although f^sc — 0-05 is found in LBGs at 
z ~ 3 (IShapley et al.ll2006tllwata et al.ll2009l) . the escape frac- 
tion of ionizing photons at z. > 4 has not been observationally 
determined yet. ^ Ono et afl (1201 0) put an upper limit on /esc at 
high redshift as /esc ^ 0.6 at z = 5.7 and /esc < 0.9 at z = 6.6. 

Fig. |3]shows the ionizing photon emissivity, dnion/dt, in the 
unit of l/s/Mpc-' comparing with various observations. We 
set /esc = 0.2 for all the observed data as well. The observed 
ionizing photons are derived from galaxy UV LFs at z = 4-7 
(lYoshi da et al. 2006; Bouwe ns_et_al. 2008; Ouchi et al. 20091 
filled symbols in Fig. O by lOuchi et al.i (2009). The conver- 
sion from LF to ioniz i ng ph oton emissivity is based on the 
Eq. 5 in lOuchi etaTI (120091) . where continuous star forma- 
tion history is assumed for all galaxies. The data are inte- 
grated down to L = 0. We als o show the ionizing photo n rate 
inferred from the Lya forest ( Bolton & H aehnelt '2007*. open 
symbols in Fig. O, which is derived by combining hydrody- 
namical simulations with measurements of the Lya opacity of 
the IGM. 

As shown later, the baseline Mitaka model does not pro- 
duce enough ionizing photons for the measured Thomson 
scattering optical depth, although it can produce sufficient 
ionizing photons to reionize the Universe at z J; 8. We 
introduce two parameters to make it consistent with vari- 
ous reionization data. One is the normalization factor X 
which renormalize the ionizing photon emissivity by a fac- 
tor of X and the other is the normalizing redshift Zc where 
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we renormalize the ionizing photon emissivity by a factor 
of X only at z > Zc from the baseline Mitaka model. In 
Fig. |3] we show models with {X,Zc) = (1.0,0.0) (baseline 
model), (50.0, 10.0), and (50.0, 10.0). Hereafter, we name 
(X,Zc) = (1.0,0.0), (50.0, 10.0), and (100.0, 10.0) model as the 
baseline model, the lower-Pop-III model, and the upper-Pop- 
III model, respectively. Fig. [3]shows the contribution of Pop- 
III stars in the baseline model. 

Although the overall behavior of our ionizing photon emis- 
sivity model at 2 < z < 7 is similar to that of the observed 
data, we overestimate the ionizing photon emissivity about a 
factor of 2 comparing to the data. We note that the observed 
data are very sensitive to the assumed spectral in dex of ioniz- 
ing emission a (see Eq. 5 in lOuchi et al.ll20()9h which is not 
well determined, a is set to be 3.0 for the data in the Fig. [3] 
Once we set a = 1.5, it will double the data derived from the 
galaxy UV LF and the Lya forest opacity. 

Fig. [3] also shows the required dNion/dt to ba lance recom- 
bination of hydrogen IGM based on the model of lMadau et al.l 
(il999i) . 



dni. 



10'*^'^C(l+z)' [s-'Mpc"^] 



dt 



frec(z) 



(2) 



where nj^ is the number density of hydrogens at present-day, 
frec(z) is the recombination time scale at z, and C = {njf)/njj is 
a time dependent volume- averaged clumping factor We show 
the cases of C = 1,3, and 10. C = 1 corresponds to the ho- 
mogeneous Universe. Each of our models predict that there 
was a sufficient ionizing photon energy budget to ionize the 
Universe at z < 7 - 8 depending on C. However, the base- 
line model can not ionize the Universe before the epoch of 
z = 8. If the Universe is reionized instantaneously, it might 
be occurred at z = 10 .6 zt 1.2 inferred from the WMAP data 
dKomatsu et al.ll20Tl1) . Even if the cosmic reionization is not 
instantaneous, we may need more ionizing photons above 
z ^ 8 t han the ionizing photon emissivity limit by .Madau et al] 
(fT999h . 

3.2. Probing the Cosmic Reionization History 

Important indicators of the reionization history are the 
Thomson electron scattering o ptical depth and t he ne utral 
fraction of hydrogens. Following lBarkana & Loebl(l2001h . we 
compute the reionization history of the Universe. The ioniza- 
tion equilibrium equation in terms of volume filling factor of 
HII region Qhii is given as 



dQ 



dt 



1 dn-ioi 
n% dt 



(3) 



where f is the cosmic time, = X//n^ is the number den- 
sity of hydrogen at present-day, is the number density of 
baryons at present-day, and Xh = 0.76 is the hydrogen mass 
fraction, dnion/dt is the ionizing photon production rate (see 
Fig. |3]l, as = 2.6 x lO"'-' cm-'s"' is the hydrogen recombina- 
tion rate, and a{t) is the scale factor. The recombination time 

scale frec(z) in Eq. |2]is given a s a(t)^ /a fiCn%. 

Equation|3]is solved as (see lBarkana & Loebll2001l for de- 
tails) 



2hii(zo) ■ 



dz 



1 dnjoj 
n% dt 



(4) 



where a constant C is assumed and dt /dz is calculated by the 
Friedmann equation in the standard flat universe cosmology 
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Fig. 4. — Thomson scattering optical depth of the Universe. We set /esc = 
0.2 and C = 3.0. Solid, dashed, and dotted curve con'esponds to the base- 
Une, the lower-Pop-III, and the upper-Pop-III model, respectively. Shaded 
region shows the 7-y ear WMAP results with l-cr errors Tc = 0.088 it 0.015 
(Komatsu et aU201 II) . 



C)=(0.2, 3.0) 

(X, zj=(1.0,0.0) — 
(50.0, 10.0) - ■ 
(100.0, 10.0) ■■■■ 



as 



dt 



dz {\+z)H,,^nM(i+z?+nf 



(5) 



Once gffli reaches 1, there is no neutral hydrogens and ioniz- 
ing photons are saved in the intergalactic space. The function 
/^(z,Zo) takes into account of recombination and given as 



2 

F(z,zo) = --^^C[/(z)-/(zo)], 



where /(z) is defined as 



(l + z)3+- 



M 



(6) 



(7) 



In this paper, we set C = 3 .0 motivated bv lPawlik et al.l (120091) . 
We do not change C as a function of redshift. We also as- 
sume QhcII = Ghii and neglecting the electron populations 
from QhcIii region to compute the optical depth of electron 
scattering. Then, the number density of free electrons at z is 
written as 



ne{z)=\Qmi{z)XH + 
\ + 3Xh 



eHeIl(l-XH)\^„^j^^^3 



Qm,{z)n'i{\+z?. 
The optical depth to the Thomson electron scattering is 



(8) 



Te(ZQ) = 



dz—aTne(z), 
dz 



(9) 



where ctj- is the Thomson cross section and dl/dz is the cos- 
mological line element given as 



dl dt 



dz dz (l+z)Hoy^nM(l+zy+nA 



(10) 



for a standard flat universe cosmology. 

Fig. m shows the Thomson scattering optical depth of the 
Universe. We also show the 7-vear WMAP data which reports 
Te = 0.088 ±0.015 dKomatsu etalj|2011,) . The baseline model 
can not reproduce the WMAP data due to the insufficient ion- 
izing photons at z > 8 (see Fig. O, although it can reionize the 
Universe at z < 8. The lower-Pop-III and upper-Pop-III model 
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Fig. 5. — Hydrogen neutral fraction of the Universe. We set /esc = 0.2 
and C = 3.0. Solid, dashed, and dotted curve corresponds to the base- 
line, the lower-Pop-III, and the upper-Pop-III model, respectively. We 
also show Gunn-Peterson trough test constraint s on the neutra l fraction 
from quasars ^Mesinger & Haiman 2007; Mortlock et alj 1201 ID . a GRB 
ITotani et aU2006l) . and LAEs jOuchi et aU2010l) . 



reproduces the lower-limit and upper-limit of the WMAP data, 
respectively. Therefore, we need 50-100 times more ionizing 
photons at z > 10 to be consistent with the WMAP Thom- 
son scattering optical depth measurement than those naturally 
expected from a semi-analytical galaxy formation model. 

Fig. |5] shows the hydrogen neutral fraction of the Uni- 
verse, Xhi = 1 - Qhti ■ Gunn-Peterson trough test constraints 
dGunn & Peterson 1965) from quasars (Mesinaer & Haiman 
|2007; Mortlock et al. 2011), a GRB (Tot ani et al . 2006), and 
LAEs (Ouchi et al. 2010) are also sho wn. We do not show the 
constraints from quasars at z < 6 (e.g. i Fan et alj|2006) . This 
is because we do not include the neutral hydrogen distribu- 
tion which is not observationally well understood. However, 
this does not affect our argument, since we trace the global 
reionization history. We note that the neutral hydrogen distri- 
bution is important to trace the ionization state of small scale 
structures. 

As shown in Fig. |5] the three models well trace the history 
of the neutral fraction of the Universe. Since we renormalize 
the ionizing photon emissivity above Zc for the lower-Pop-III 
model and the upper-Pop-III model, ionizing photon emissiv- 
ity suddenly decreases at z = Zc and the Universe starts to re- 
combine ionized hydrogens again. This is because these two 
models show a bump feature at z ^ 7. 

As shown in Figs. |4]and|5] 50-100 times more ionizing 
photons are required at z > 10 to be consistent with the Thom- 
son scattering optical depth measurement comparing to the 
baseline model which successfully reproduces the various ob- 
servational data of galaxy formation and generates sufficient 
ionizing photons to reionize the Universe at z < 8. Observa- 
tionally, it is also known that the estimates of the Thomson 
optical depths and contribution of galaxies to reionization are 
higher than those in observational studies of galaxies detected 
atz<8(e.g . Stark et al. 2007; Charv 2008; Oesch et al. 2009; 
Ouchi et al., |200^ PawHk et al. 2009; Bu nker et al.i i2010 ; 
Labbe et alj boiot [Robertson 2010; Bouw ens et all 12011a , 



20121) . The enhancement of ionizing photon emissivity might 



be required at high redshift to be consistent with the WMAP 
measurement of the Thomson scattering optical depth. 

There are various possible ingredients which play an im- 
portant role in generating sufficient ionizing photons for the 
WMAP data. These are the steeping of the faint-end slope in- 
dex of LF (.Bouwens et al...201Z) . a smaller clumping factor 



dBolton & Ha ehnem 2007HPawhk et alj|2009h . a larger escape 
fraction (jYaiima et alj 201 Ih. a harder initial mass function 
at higher redshift (McKee & Tan"2008'), a contribution from 
early Pop-Ill stars (Cen 2003), and feedback from accreting 
black holes (Ricotti & Ostriker 2004; Mirabel et al. 201 1). 

To study Pop-Ill stars' contribution to the NIR background, 
we assume that required ionizing photon comes from Pop-Ill 
stars. This assumption will maximize the Pop-Ill stars' con- 
tribution to the NIR background. Therefore, we renormalize 
the Pop-Ill emissivity by a factor of X at z > Ze- As shown 
later, this does not significantly affect our EBL and gamma- 
ray optical depth models. 

4. EXTRAGALACTIC BACKGROUND LIGHT 

We compute the background intensity I{i^i),zq) at redshift 
Zo and frequency i/q which is the integrated radiation from 
all sources between z = Zo and the maximum redshift of the 
source distribution, Zmax- Their contr ibution to the lo cal back- 
ground intensity is given by (see e. g. iPeacockll 1 9991) 



2 /'-i-max 
A-TT 



dz—j(v,z), 
dz 



(11) 



where j(v,z) is the galaxy comoving volume emissivity at 
redshift z and frequency i' = i^q(1+z)- We set Zmax = 20. The 
galaxy comoving volume emissivity is calculated by combin- 
ing our CSFH and stellar population synthesis models. 

For the units in the proper volume, from Eq. [TT] the specific 
radiation energy density (in the unit of erg/s/cm^/Hz) in the 
proper volume is 



4-7r 

pi'^o,zo) = — (1 + zo)^I(i^o,zq) 
c 



(12) 



where a factor of (1 +zo)^ corrects the energy density from 
local to Z(). The photon proper number density is 



dn(eo,zo) p(va,Zo) 



deo 



(13) 



where eo = hpiyQ is the photon energy and hp is the Planck 
constant. 

Fig. |6] shows plots of the proper photon number density 
based on the baseline, upper-Pop-III, and lower-Pop-III mod- 
els. When we increase the Pop-Ill population, we increase 
their stellar radiation only. Thus, there is a discontinuity at 
^ 0.2 eV (^ 6.2 fim). Proper photon number density in- 
creases from high redshift to z ^ 1 - 3 where the CSFH has 
a peak and decreases toward the local Universe. Once we in- 
crease the Pop-Ill component, the UV photon density at high 
redshift becomes comparable to that at z = 1 . 

The local background intensity reported in Fig. |7]as a func- 
tion of wavelength A is calculated from Equation[TT|by setting 
Zo = 0. We show the baseline model here . We also show other 
theoretical models for the comparison (iKneiskeet all 12001 



Stecker et al] 120061; [F ranceschini e t"aLl 120081; i Gilmore et al 
2009 tlFinke et al.ll2010^,Kneiske & Dolell20TalGilmore et al 
20121) . The current measurements of the EBL and the 
integrated brightness of galaxies are also shown. Detailed 
predictions for the proper photon number density and the local 
background intensity is publicly distributed on the Web site 
http : / / www . s lac . Stanford . edu/~yinoue/ Dow nload. 

Overall shape of our EBL model is consistent with the ob- 
servational data. Our model is in a good agreement with the 
observations by Pioneer 10/11 (.Matsuoka et al...201 L open 
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Fig. 6. — Proper volume photon number densities multiplied by the photon energy e as a function of redshift. Top-left, top-right, and bottom-left panel 
conesponds to the baseline model, the upper-Pop-III model, and the lower-Pop-III model, respectively. For the upper-Pop-III and lower-Pop-III models, we only 
renormaUze their stellar contribution. SoUd, dashed, dotted, dot-dashed, double dot-dashed, and triple dot-dashed curve corresponds to the proper photon density 
at z = 0.1, 1.0, 3.0, 6.0, 8.0, and 10.0, respectively. 



pentagon symbols in Fig. |7]l which directly measured the 
EBL from the outside of the zodiacal region. Our model also 
does not violate th e limit from the gamma-ray observation 
(lAlbert et alJl2008l) . 

Comparing to other models, we tends to predict more pho- 
tons at A < 0.4/im and less photons at A > 0.4/im. At 
UV range, other models except for Steckeretal. (2006) 
are consistent with the GALEX data (Xu et al. 2005), while 
ours are consistent with the HST data ( Gardner et al. 2000). 
Both observational data points show the integrated galaxy 
counts down to zero luminosity. The reason why we pre- 
dict more U V photons is due to the treatment of the dust 
obscuration. iSomerville et alj (120121) sho wed that a model 
with the Calzetti law (^Calzettietal.ll2000^ ■ which we adopt, 
predicts more UV photons comparing to models with the 
multi-dust component model which Gilmore et al. (20l3); 
ISomerv ille et al. (2012) a dopted. W e also note that the data 
points bv .Madau & Pozzettir ( l200(il) are the integrated galaxy 
counts down to the dete ction Umit of the HST. Therefore, 
iMadau & Pozzettil (l200^ gives the lower-limit of the EBL. 

At 0.4/im < A < 10 /im, ou r model is reasonably in a good 
agreement with iKneiske & Dold (12010 ) which predict the 
lower-limit of the EBL. The reason why we predict less pho- 
tons at A > 0.4/im comparing to other models is that our 
CSFH is a factor of 3 lower than that of iHopkins & BeacomI 
12006) which are used in most of previous studies. 

Fig. [8] shows the Pop-Ill EBL from the baseline, upper- 
Pop-III, and lower-Pop-III models. Each model shows the 
Pop-Ill contribution is < 0.03 nW/m^/sr which is less than 



0.5% of the total NIR background radiation. It is difficult to 
explain the IRTS data (Matsumoto et al. 2005) by the Pop-Ill 
population even with the highest level of the arrowed Pop- 
III contribution from the reionization constraints. Moreover, 
even if the trapped ionizing photons in a galaxy is converted 
to Ly-a photons, the NIR flux from Pop-Ill stars will increase 
only 15% in the case B recombination. Therefore, it is hard to 
extract the Pop-Ill star information from the NIR background 
flux measurements. Since we renormalize the Pop-Ill popula- 
tion at z > Zc = 10, two peaks appear in the upper-Pop-III EBL 
and the lower-Pop-III EBL. 



5. GAMMA-RAY ATTENUATION 

5.1. Gamma-ray Opacity 

The redshift dependent background photon number density 
is now obtained by following ^ We are able to compute 
the gamma-ray opacity for pair production interaction s. The 
cross section of pair production process is given as (iHeitleii 
IT954h 



a^^(£^,e,0)=^(l-/?2) 



2/3(/32-2) + (3-/3'*)ln 



1+/ 
1-/ 



, (14) 



where e is the energy of the background photon, is the 
energy of colUding high energy photon, 9 is the angle between 
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Fig. 7. — The EBL model for the baseline model is shown by solid curve. For comparison, the EBL models bv |Kneiske et alj <2004t thin dashed), 'Stecker et al. 
(2006, dotted), Franceschini et al. 1 2008, dot-dashed), Gilmore et al. 1 2009, thin double dot-dashed), Finke et al. (2010, triple dot-dashed), Kneiske & Dole 1 201^ 
thick dashed), and Gilmore et al. (2012, thick double dot-dashed) are shown as indicated in the figure. The integrated brightness of galaxies (minimum EBL; 
filled symbols) and current measurements of the EBL (open symbols) are shown as indicated in the figure. References for the integrated brightness of galaxies 
are HST {Madaa & Pozzetti 2000; Gardner et al. 2000), ISO (Elbaz et al. 2002), Spitzer fPapovich et al. 2004; Fazio et al. 2004; Dole et al.''2006; Frayer et 
2006), GALEX l Xu et al. 2005). and Henhel I Berta et al. 201 1 ). References for the current EBL measurements are DIRBE (Wrisht & Reese 2000; Wright 20O 
Levenson et al. 2007; Levenson & Wrisht 2008; Cambresv et al. 2001; Dwek & Arendt 1998; Goriian et al. 2000; Finkbeiner et al. 2000; Hauser et al. 199: 
Lagache et alji200Q), HST »Bernsteini2007i ; iBrown et al.ii2000) . Pioneer I Mats uoka et al .1,2011.) , IRTS I Matsumoto et alji2005i) . AKARI I Matsuur a et al.,i20Hl 
and Voyser 4Edelstein et alJ200(D . The upper limit from TeV gamma-ray observation by MAGIC j Albert et al. 2008} is also shown by solid line with arrows. 



r IcL 


i\ 




t| 1 ; 














All:(X,Ji,H1.0, 0.0) 

Pop-Ill: (X,z^)=( 1.0, 0.0) 

Pop-Ill: (X, z^)=(50.0, 1 0.0) : 

Pop-Ill: (X, z^)=(I(X).0, 10.0) 











10 
Xlflml 



Fig. 8. — Same as Fig. [7] but showing the Pop-III contribution to the EBL. 
Thick solid curve shows the total EBL. Thin-solid, dashed, and dotted curve 
corresponds to the baseline, the upper-Pop-III and the lower-Pop-III model, 
respectively. 



the colliding photons, and /3 is 

/?= Wl- ^ ^f^""" ■ n = cos9. (15) 



Then, the gamma-ray opacity is maximized for photon energy 
is as shown in Eq. [T] 

The optical depth for a high energy photon E^^ propagating 
the universe from a source at redshift Zs to an observer at zo is 



dz I 



de 



dl I— fi 
^'dz~l~ 



x^^^a.,^(E.,(l+zXe,0), (16) 

where eth is the threshold energy for pair production process 
and given as 



eth ■■ 



£^(l+z)(l-/x) 



(17) 



££^(1 -cos6') 



Fig. |9] shows the optical depth for 77 interaction as a func- 
tion of the observed gamma-ray energy E~, from z = 0- 10 
comparing with va r ious previous models (Kneiske et al.l 2004t 
Franceschini et al. 2008^, iFinke et all 120 lOt iGilmore et alJ 
2012; Inoue et al. 2010a). We calc ulate it usi ng Eq. [T6lby 
setting Zo = 0, while the model by llnoue et al. (2010a) has 
set Zo = 4. We have not considered the absorption by the 
CMB photons which is negligible in the parameter range of 
our interest at 1 GeV - 100 TeV. We show here the baseline 
model only, since the Pop-III contribution is negligible (see 
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Fig. 9. — Gamma-ray optical depth for the observed gamma-ray energy from z = 0.1, 0.3, 1.0, 3.0, 5.0, 7.0, and 10.0. Solid curve shows the baseline model 
which is comparable to other models in this paper. Dashed, dotted, dot-dashed, double dot-dashed, triple dot-dashed curve shows the model by Kneiske et al] 
j[2004) Franceschini et al. (2008), Finke et al. (2010), Gilmore et al. (2012), and Inoue et al. (2010a), respectively. Shaded region represents the measurement at 
the 95% confidence level of the gamma-ray opacity by Fermi at z Ri 1 <The Fentii-LAT Col laboration 2012). Horizontal thin dotted hne shows r-y^^ = 1. Since 
the difference among models with various Pop-Ill contents is small, we do not plot other models in this paper. 



By setting Zc 



below and §. IH. Detailed predictions for the 77 optical depth 

based on our model is pubUcly distributed on the Web site ^ 

|http : / /www . slac . Stanford . edu/~yinoue/Downl oa(±eii taLf Scient UV photons at z < 10. Moreover. Ilnoue et alj 



10, our model predi cts few UV photon s be- 
low z= 10 to keep reionization, while Ilnoue et al.l ( l2010al) has 



Our gamma-ray opacity model is consistent with the 
gamma-ray opacity m easurement by Fermi si z~ 1 at th e 
95% confidence level (iThe Fermi-LAT Collaboration! l2012h . 
Our model is also gener ally consistent wi th previou s mod- 
els (Kneiske et al. 2004; Franceschini et alj |2b08; Fink e et al] 
lIOlQciGilmore et al.n2012h at < 400/(1 +z) GeV at z < 5. 
However, opacity above Ej ^ 400/(1 +z) GeV is a factor of 
^ 2 transparent to previous models. We predict low EBL in- 
tensity at A > 0.4 /im (see Fig. |7]i where the gamma-ray opac- 
ity is maximized for the ^ 300 GeV photons (see Eq. [U. 

As we described before, redshift evolution of our MIR-FIR 
emissivity model is uncertain. EBL with longer than MIR 
is important for the absorption of gamma-rays above several 
tens TeV following Eq. [T] Therefore, this uncertainty would 
not affect to our results below 10 TeV. 

We predict higher energy cut-off in the gamma-ray spec- 
trum at ^20 GeV for sources at z > 6 comparing to that in 
Ilnoue et aTl (l2010al) which predicted ^12 GeV at z ^ 5 and 
~ 6-8 GeV at z > 8 - 10. This is because the UV radiation 
from Pop-III stars are not sufficient enough to significantly 
absorb gamma-rays at high redshift even in the case of the 
upper-Pop-111 model as discussed below. Since path length at 
z > 6 is not so long as that at lower redshift, the high emis- 
sivity due to Pop-Ill stars does not affect the opacity signifi- 
cantly. 



(12010a) did not include stars with metallicity higher than 
Z = 0.02Zq. These might be the main reason why our predic- 
tions are different. We note that the EBL contribution from 
quasars which is not included in our study is not significant as 
shown in llnoue et alj ( l2010a[) . 

Based on a gamma-ray luminosity function of blazars, 
Ferm i is expected to detect a blazars at z ^ 6 (Inoue et al] 
1201 Ih . Moreover, GRBs are also promising a s they are 
known to occur at z > 6 dKawai et al.l I2006|; [Greiner et al] 
I2009h . at least up to z 8.2 (^T anvir et all 120091) . and prob- 
ably ou t to the epoch of first star formation in the Universe 
dBrom m & Loeb 2006). For example, if a GRB 080916C-like 
event (Abdo et al. 2009) is occurred at z ^ 6 which was the 
brightest GRB event in the gamma- ray sky , CTA can clearly 
detect its gamma-ray spectrum (Ino ue et aU lin preparation!)- 

Fig. [To] shows the 77 opacity from all stars. Population I 
(Pop-1) stars (10-2-5 < Z), Pop-11 stars (lO""* < Z < lO'^^), 
and Pop-Ill stars (Z < 10""*). Since the metallicity transi- 
tion between Pop-11 stars (halo stars) and Pop-1 stars (disk 
stars') is [Fe/H]^-1 (jWhe eler et alJ[T98l lMcWimaml ll997l: 
iProchaska et al.l 120001) corresponding to Z ~ 10"^-^- 10"2 \ 
we take Z = lO"^'^ for the transition metallicity between these 
populations in this paper. We show here the baseline model. 
Although it is difficult to capture the direct signature of Pop- 
III stars in the gamma-ray spectra of high redshift sources, the 
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Fig. 10. — Same as Fig.|9] but shows the gamma-ray optical depth from each stellar populations. Solid, dashed, dotted, and dot-dashed curve shows contribution 
of all stars, Pop-1 stars, Pop-11 stars, and Pop-Ill stars, respectively, in the baseline model. Horizontal thin dotted line shows t-^.^ = 1- Since the contribution of 
Pop-Ill stars is weak, it does not appear in the panels at z = 0. 1 and 0.3. 



attenuation due to Population II stars should be observable in 
future high-z gamma-ray observation and will provide a valu- 
able probe of the evolving UV EBL in the cosmic reionization 
epoch. Moreover, once high redshift gamma-ray events are ob- 
served, we can at least put an upper limit on the Pop-Ill stars 
contribution to the EBL and may be able to distinguish the 
reio nizaiton history of th e Universe by comparing our model 
with llnoue etal] ( l2010aD . 

Fig. [TT| shows the differential of gamma-ray optical depth 
dr^^/dz with respect to redshift for the baseline, lower- Pop- 
Ill, and upper-Pop-III models. Since the Pop-Ill component 
is added at z > 10, the emissivity contribution from Pop-Ill 
is insignificant at z < 6. At z > 6, the Pop-Ill contribution to 
dr^-y/dz depends on the amount of the Pop-Ill stars. There- 
fore, there is a difference in gamma-ray opacity among mod- 
els at high redshift. This results in -3%, -10% and -20% 
enhancement of the gamma-ray opacity at 20 GeV compar- 
ing to the baseline model at z = 6, 8, and 10, respectively. At 
higher energy, the enhancement decreases. If it is possible to 
identify — 10% difference in flux, we are able to distinguish 
our models. 

5.2. Comparison with Current GeV & TeV data 

Gamma-ray astronomy has now been firmly established 
by observations of Fermi, H.E.S.S., MAGIC, and VERI- 
TAS. Further progress is anticipated in the near future by 
CTA. The most distant gamma-ray objec t detected in GeV 
and TeV band is a GRB at z = 4.35 dAbdo et all l2009h 
and a blazar at z = 0.536 dAlbert et alj l2008h . If the 
gamma-rays observed from the direction of these sources 



are created at the source redshift and are generated by 
a single component such as shock accelerated particles, 
the gamma-ray opacity of the Universe h a s already been 
constrained (see e.g. M azin & Raud I2007t iRaue & MazinI 
120081) . CTA is ex pected t o detect > 100 blazars up to 
z - 2.5 dlnoue et al.ll2 010b: In oue & CTA Consortiumll201 It 
ISol et alJIin preparatio n). Future high redshift blazar observa- 
tion will be a key to understanding the evolution of the opt.- 
NIR EBL in detail. Precise low redshift TeV sources obser- 
vation will be a key to understanding the FIR EBL, since the 
opacity is huge at high redshift. Although starburst galax- 
ies are suggested as possi ble targets for the FIR EBL study 
dPwek & Krennrichll2012l) . the internal gamma-ray absorp- 
tion will not allow photons above 10 TeV escaping from a 
galaxy (e.g. Inoue 201 lb). 

A plot of the energy at which the Universe becomes opti- 
cally thick to gamma-rays, defined where t^^ = 1 is shown 
in Fig. [T2I which is known as the Fazio-Stecker relation 
dFazio & Steckei^'igVO*). We show here our baseline model 
and other models (Kneiske et al. 2004; Franceschini 1 



2008t iFinkeet all 12010: IGilmore et alJ 120 12t ilnoue 



et aD 
etall 



2010al) . Since other models do not compute the opacity at 
z = - 10, we do not show the opacity of each model at 
the outside of the redshift range of each paper We also 
plot the maximum photo n energy bins of several blazars (see 
'Finke & Razzaaud'2009', for a list and references) and GRB 
080916C(Abdo et al. 2009). 

Most of blazars are highly attenuated by the EBL, since 
several sources are considerably above the t^^ for all models. 
The highest energy photons from the GRB 0809 16C constrain 
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Fig. 1 1 . — Differential of gamma-ray optical depth dr^-y/rfz for the observed gamma-ray energy E-^ with respect to redshift i from z = 1.0, 4.0, 6.0, 8.0, and 
10.0. Solid, dashed, and dotted curve shows the basehne, lower-Pop-III, and upper-Pop-III model, respectively. Each curve overlaps each other at z = 1 and 4. 
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Fig. 12. — Horizon of the Gamma-ray Opacity where t^^i = 1, known as the 
Fazio-Stecker relation (Pazio & Stecker 1970). The baseline model is shown 
by solid curve. Dashed, dotted, dot-dashed, double dot-dashed, triple dot- 
dashed curve show s the model bv[Kn eiske et al. (2004) Franceschini et al. 
f200^. 'Finke et al.' l'20T(T). 'Gil moreet al.i >2012i) . and jinoue et al. (2010a), 
respectively. Thin dashed curve show the case of r^-y = 5 for the baseline 
model. The filled data points are the observed maximum gamma-ray photon 
energy bins from blazars (Finke & Razzaque 2009) and the GRB 080916 C 
(Abdo et alj|2009i) . Since other papers do not cover the opacity at z = - 10, 
we do not show the opacity of each model at the outside of the redshift range 
of each paper. 



the EBL at high redshifts. Our model predicts that the Uni- 
verse is transparent below 20 GeV even at z > 4. 

Figs. [T3] and [T4l show the observed spectra of TeV blazars 
at z < 0.15 and z > 0.15, respectively, and their deabsorbed 
spectra. If TeV emission of these sources are originated from 
shock accelerated electrons or protons, the intrinsic photon 



index Tint = 1.5 is expect ed to be th e hardest in the diffusive 
shock acceleration (Malk ov & O'C Drury 2001). Although 
the results are generally consistent with a minimum intrin- 
sic Tint = 1.5, some sources such as lES 0229-20 and lES 
1 101-232 shows harder spectra. Moreover, it seems that there 
are two components in some sources, since spectra starts to 
harden from several hundreds GeV. 

Some authors recently claim the existence of the secondary 
cascade component in these energy band from very high 
energy cosmic-rays or gamma-rays ( e.g. Essey et al. 201lt 
Essey & Kusenkoll2012l : iMurase et al.ll2012 : Ahai'onian et al] 
20121) . These components will be a key to understanding the 
IGMFs. Once future CTA observations bring us spectra of 
these sources with a high energy resolution, we will be able 
to search such signature by comparing with the deabsorbed 
blazar spectra. 

The EBL absorption signature is not seen in the spectrum 
of the extragalactic gam ma-ray background (EGB) spectrum 
jAckermann et alj|201 II) . although it is naturally expected i f 
the origin is cosmological ( Ilnouel201 latllnoue & Iokall2012h . 
Inoue & loka (2012) has recently showed that the EGB mea- 
surement at > 100 GeV (Ackermann et al. 201 1) violates the 
EGB itself at < 100 GeV (.Abdo et al...2010.) by considering 
known sources' contribution. The upper hmit on EGB is de- 
rived from the cascade emission from EGB above 100 GeV 
not to exceed the current EGB measurement. Even if our 
EBL model is adopted, the violation remains. Detailed stud- 
ies of gamma-ray spectra of extragalactic sources are required 
to probe this violation. 

6. CONCLUSIONS 
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We have developed models for the EEL over the redshift 
range z = 10toz = 0on the basis of a semi-analytical model of 
hierarchical galaxy formation, into which Pop-Ill stars were 
incorporated in a simplified fashion. Our baseline model is 
consistent with a wide variety o f observational data fo r galax- 
ies below z ^ 6 (iNagashima & Yoshii 2004; Kobayashi et al.l 
I2OO7I lIoToh . and is also capable of reionizing the Universe 
by z < 8. However, in order to account for the Thomson 
scattering optical depth measured by WMAP, the ionizing 
photon emissivity is required to be 50-100 times higher at 
z > 10. This is line with recent observations of galaxy can- 
didates at z ^ 8, as long as the contribution from faint galax- 
ies below the se nsitivity of current telescopes is not large (e.g. 
iBouwens et al.l l2012). The "missing" ionizing photons may 
possibly be supplied by Pop-Ill stars forming predominantly 
at these epochs in sufficiently small galaxies. 

The EBL intensity at z = in our model is generally not 
far above the lower limits derived from galaxy counts. Our 
model is also in good agreement with the data from Pioneer 
dMatsuoka et al.ll20Tl ) directly measured from outside the zo- 
diacal region. The Pop-Ill contribution to the NIR EBL is 
< 0.03 nW/m^/sr, less than 0.5 % of the total in this band, 
even at the maximum level compatible w ith WMAP mea- 
surem ents. The putative NIR EBL excess dMatsumoto et al.l 
I2005h . which also conflicts with the upper limits from gamma- 



ray observations (lAharonian et al.ll2006ah . may have a zodia- 
cal origin rather than Pop-Ill stars. 

Up to z ^ 3-5, the 77 opacity in our model is compa- 
rable to that in the majority of previously published mod- 
els (Kneiske et al. 2004; Franc eschini et al.H2008l: iFinke et akl 
I2OIO; Gilmor e et al. 2012) below ~ 400/(1 +z) GeV, while 
it is a factor of 2 lower above this energy. The Universe 
is predicted to be largely transparent below 20 GeV even at 
z > 4. 

Estimates based on the observed gamma-ray luminosity 
function of blazars show that Fermi may detect blazars up to 
z - 6 (Inoueetal.l l20TTl). CTA may possibly d etect GRBs up 
to similar redshifts (llnoue et aLllin preparationh . However, the 
contribution of Pop-Ill stars may be difficult to discern in the 
attenuated spectra of high-redshift gamma-ray sources, even 
at the highest levels allowed by the WMAP constraints. Nev- 
ertheless, the signature of Population II stars is expected to be 
observable in high-z gamma-ray sources, providing a unique 
and valuable probe of the evolving EBL in the rest-frame UV. 

We thank Floyd Stecker and Alberto Dominguez for help- 
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the Research Fellowship of the Japan Society for the Promo- 
tion of Science (JSPS). SI is supported by Grants-in-Aid Nos. 
22540278 and 24340048 from MEXT of Japan. 
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